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ABSTRACT 

Context. The prediction of the spins of the compact remnants is a fundamental goal of the theory of stellar evolution. 
Aims. Here, we confront the predictions for white dwarf spins from evolutionary models including rotation with observational con- 
straints. 

Methods. We perform stellar evolution calculations for stars in the mass range 1. . . 3 M Q , including the physics of rotation, from the 
zero age main sequence into the TP-AGB stage. We calculate two sets of model sequences, with and without inclusion of magnetic 
fields. From the final computed models of each sequence, we deduce the angular momenta and rotational velocities of the emerging 
white dwarfs. 

Results. While models including magnetic torques predict white dwarf rotational velocities between 2 and 10kms~', those from the 
non-magnetic sequences are found to be one to two orders of magnitude larger, well above empirical upper limits. 
Conclusions. We find the situation analogous to that in the neutron star progenitor mass range, and conclude that magnetic torques 
may be required in order to understand the slow rotation of compact stellar remnants in general. 
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1. Introduction 

During the last decade, rotation and rotationally induced trans- 
port processes have been a major focus of the theory of mas- 
| sive star evolution (Maeder & Meynet 2000, Heger et al. 2000). 
Evolutionary models without magnetic field induced internal an- 
gular momentum transport predict spins of newly born neutron 
stars which are 1...2 orders of magnitude above those deduced 
from the youngest Galactic pulsars (Ott et al. 2006). However, 
the magnetic torques proposed by Spruit (1998, 2002) pro- 
vide enough coupling between the rapidly rotating core and the 
slowly rotating envelope in post-main sequence stars to produce 
the right amount of core spin-down (Heger et al. 2005, Petrovic 
et al. 2005). 

In models of low mass stars, the role of non-magnetic ro- 
tational transport processes for the determination of the white 
dwarf spin is rather negligible, which leads to rapidly rotating 
CO-cores (Langer et al. 1999, Palacios et al. 2003, 2006; cf. be- 
low). The huge shear above these rotating CO-core during the 
thermally pulsing AGB phase inhibits the s-process by swamp- 
ing the 13 C-pocket with the neutron poison 14 N (Herwig et al. 
2003, Siess et al. 2004). Analogous to the situation in massive 
stars, it appears that such rapidly spinning cores in low mass 
star models also contradict direct observations of compact stel- 
lar remnants. In a recent analysis of the Ca line profiles of a large 
sample of DA white dwarfs, Berger et al. (2005) concluded that 
their rotational velocities are generally below 1 km s~ 1 , which is 
the smallest upper limit derived from spectroscopy so far. They 
conclude that the predicted white dwarf spin from the rotating, 
non-magnetic evolutionary models of Langer et al. (1999) can 
be ruled out. This conclusion is confirmed by Kawaler (2004), 
who argues that the rotation rate of ZZ Ceti pulsators is even at 



least one order of magnitude slower than Berger et al.'s upper 
limit. 

Internal magnetic torques have been suggested as an agent 
to spin down the cores of white dwarf progenitors during the 
giant stage (Spruit 1998). In this paper, we investigate whether 
magnetic torques, as computed in Spruit (2002), are able to alle- 
viate the problem of slowly spinning white dwarfs. Eggenberger 
et al. (2005) already used the same physics to compute the an- 
gular momentum evolution of a solar mass star, with the result 
that the flat internal angular velocity profile of our Sun could be 
recovered. 

Magnetic torques may not be the only mechanism to pro- 
vide a spin-down of the stellar core. E.g., Charbonnel & Talon 
(2005) invoked angular momentum transport though gravity 
waves (Zahn et al. 1997) to explain the slow rotation of the 
Solar core. However, here we concentrate on the magnetic torque 
mechanism, as it gives promising results for massive stars and 
can be readily applied to the low mass regime. 



2. Method 

We use a 1-D hydrodynamic stellar evolution code (Yoon et al., 
2006, and references therein). It includes diffusive mixing due to 
convection, semiconvection (Langer et al. 1985) and thermoha- 
line mixing as in Wellstein et al. (2001). This code includes the 
effect of the centrifugal force on the stellar structure, and time- 
dependent chemical mixing and transport of angular momentum 
due to rotationally induced instabilities (Heger et al. 2000). We 
also include chemical mixing and transport of angular momen- 
tum due to magnetic fields (Spruit 2002) as in Heger et al. (2005) 
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and Petrovic et al. (2005). We assume rigid rotation in convec- 
tion zones. 

We calculate the evolution of solar metallicity stars start- 
ing rigidly rotating at the zero-age main sequence, for initial 
masses of 1.0 M , 1.5 M , 2.0 M and 3.0 M , with and without 
the effects of magnetic torques. The initial equatorial velocities 
of these models were chosen to be 2, 45, 140 and 250 km s -1 
(Tassoul, 2000). For the solar mass models, we assume that the 
angular momentum loss due to magnetic braking connected with 
the solar-type wind is already over once we start our calcula- 
tions (cf. Eggenberger et al. 2005). All but our 3 M Q models 
evolve through the core helium flash. Throughout the evolution 
of all models, the mass loss rate of Reimers (1975) was used. 
The choice of the AGB mass loss rates does not affect our re- 
sults significantly. The calculations were stopped after at least 
five and at most 28 thermal pulses, beyond which the CO-core 
angular momentum is not expected to change significantly (cf., 
Sect. 3). 

We compute the surface rotational velocities of the white 
dwarfs emerging from our model sequences as follows. First, we 
fix the final white dwarf mass according to the initial-final mass 
relation from Weidemann (2000). These white dwarf masses 
(cf . Table 1 ) are generally slightly larger than the final CO-core 
masses of our models. However, as we shall see below, the total 
angular momentum of the CO-core does not change any more 
during its growth on the TP-AGB, and therefore, even though 
we correct our final CO-core masses in the mentioned way, our 
final angular momenta do not need to be adjusted. 

The surface rotational velocities are then calculated assum- 
ing rigid rotation inside the white dwarf. This is well justified 
due to the short timescale of angular momentum redistribution 
(~ 10 8 yr) inside even the non-magnetic models (Yoon & Langer 
2004). We then use the gyration radii, H, for polytropes de- 
rived by Motz (1952). For a white dwarf with a polytropic index 
n = 1.5, mass Mwd, radius Rwd (we employ the white dwarf 
mass-radius relation of Hamada & Salpeter, 1961), and angular 
momentum / tot , it is (Table 1 in Motz, 1952): 



H 



0.20502 



(1) 



With <d = J/I and /wd = H 2 M^d we obtain the surface velocity 



V rot ,WD = ^WD*WD = — *WD = a2()5()2/;wDMwD > 



(2) 



where the angular velocity of the white dwarf. We have 

tested this procedure by computing rotating cold white dwarf 
models with our stellar evolution code and find it accurate to 
within 10% for the slowly rotating cases considered here. 

3. Results 

The evolutionary tracks of our four magnetic models in the HR- 
diagram are shown in Fig.fT] while Fig. [2] displays their evolu- 
tion in the T c - p c -plane. The tracks of our non-magnetic se- 
quences would be indistinguishable on these plots. Fig.|2]shows 
that during the late AGB evolution, the CO-cores of our models 
are degenerate and cooling. Their central densities are close to 
that of cold white dwarfs in the the considered mass range, i.e. 
p c 3. ..6 10 6 gcnr 3 forMwo = 0.55 M ... 0.7 M Q . 

In order to visualize the angular momentum distribution in 
our models, we show diagrams using the integrated angular mo- 

mentum J(M r ) = J () j(M r )dM r divided by M/ , as a function 
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Fig. 1. Evolutionary tracks of the four sequences including mag- 
netic torques in the Hertzsprung-Russell diagram (cf. Table 1). 
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Fig. 2. Evolution of the four sequences including magnetic 
torques in the central density (p c ) versus central temperature (T c ) 
plane 



of the mass coordinate M r . While a homogeneous, rigidly ro- 
tating body results in a horizontal line in such a diagram, the 
/(M,)-profiles from different evolutionary stages trace the flow 
of angular momentum through the mass shells, since J{M r ) and 
J{M r )M r remain constant in a given mass shell if no angular 
momentum is transported through this shell. 

Fig. [3] displays /-profiles from six evolutionary stages for 
our 1.5 M Q model, both with and without the effects of mag- 
netic torques. For the non-magnetic model, the deepest extension 
of the convective envelope is down to about 0.2 M . Therefore, 
the helium core is built-up starting with the initial angular mo- 
mentum of the inner 0.2 M , i.e. about 10 48 cm 2 s . During the 
growth of the helium core on the RGB, material from the for- 
mer bottom of the convective envelope is incorporated, which 
is very slowly rotating and which, in its inner parts, has a spe- 
cific angular momentum several orders of magnitude below that 
of the material in the helium core. Thus, effectively, the helium 
core grows in mass without gaining angular momentum. This 
feature is continued for the CO-core during the AGB-evolution. 
As a consequence, the /-profiles in Fig. [3] follow the line of 
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10 48 cm 2 s- 1 



in the mass range 0.2 M < M r < 0.6 M 
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Table 1. Initial mass Mj, initial equatorial velocity u rot> j, adopted 
white dwarf mass Mwd and radius Swd final core angular mo- 
mentum / Cj f, final mass averaged specific core angular momen- 
tum j c j and expected equatorial rotation velocity (iw.wd) from 
Eq. 2, for our magnetic (upper 4 lines) and non-magnetic mod- 



0.5 : 



els. 


Mi 








icf/10 46 


Jcf/10 13 


frot.WD 


M 


kms 1 


M 


10- 2 Ro 


gem 2 s" 1 


cm 2 s 1 


kms" 1 


1.0 


2 


0.550 


1.25 


5.00 


4.54 


2.6 


1.5 


45 


0.575 


1.20 


5.83 


5.07 


3.0 


2.0 


140 


0.600 


1.19 


8.05 


6.71 


4.5 


3.0 


250 


0.680 


1.08 


15.4 


11.3 


9.1 


1.0 


2 


0.550 


1.25 


8.07 


7.33 


4.2 


1.5 


45 


0.575 


1.20 


110 


95.6 


56 


2.0 


140 


0.600 


1.19 


329 


274 


180 


3.0 


250 


0.680 


1.08 


458 


336 


220 



during the post main sequence evolution. Other processes than 
convection are clearly negligible in this case. The specific angu- 
lar momentum in the inner 0.6 M core of this model has been 
lowered by about a factor 6, from the ZAMS to the white dwarf 
stage. For the other non-magnetic models, this factor is rather 
similar. 

For the magnetic 1.5 M sequence, Fig. [3] shows a drop of 
core angular momentum already during core hydrogen burning. 
This is a consequence of the magnetic fields enforcing close-to- 
rigid rotation. However, the main drain of core angular momen- 
tum occurs between core hydrogen exhaustion and core helium 
ignition, i.e. during the RGB phase. The total loss of core an- 
gular momentum here is about a factor of 20 larger than in the 
non-magnetic case, i.e. a factor of 120 in total. Fig. [4] which 
compares the final angular momentum distributions of our mag- 
netic and non-magnetic models, shows that the situation is sim- 
ilar for all the studied cases, except for the solar mass models, 
where the initial angular momentum was already rather low. 



4. Discussion 

Fig.[5]summarizes our results by showing the initial and the final 
average specific core angular momentum of our model. It also 
contains the corresponding result form Heger et al. (2005) for 
a 15M Q star. Observational constrains from white dwarf spec- 
troscopy, pulsation analyses of ZZ Ceti stars, magnetic white 
dwarfs, and young neutron stars are also plotted on this figure. 

In the initial mass range considered here, as well as in the 
massive star regime, it is clear that the non-magnetic models 
fail by a large margin to comply with the observations. Only 
for solar mass stars, where magnetic braking due to the solar- 
type wind leads to a low initial angular momentum, could the 
non-magnetic model be marginally consistent with the empir- 
ical data. However, our assumption was that the slow surface 
rotation of these stars early during hydrogen burning has slowed 
down the rotation of the whole stellar interior. Eggenberger et al. 
(2005) showed that exactly this might have been performed by 
the B-fields produced due to the Spruit (2002) dynamo. 

The magnetic low mass models all fall below the spectro- 
scopic limit on white dwarf rotation of Berger et al. (2005). 
Towards more massive stars, the agreement with observations 
seems good. The neutron star spins are well recovered, and the 
range of spins of the magnetic white dwarfs (which are presumed 
to have more massive progenitors than typical white dwarfs; 
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Fig. 3. Integrated angular momentum J{M r ) = J Q j(m)dm di- 
vided by M, 5/3 as a function of the mass coordinate, at different 
times for the magnetic (full drawn red lines) and non-magnetic 
(dashed blue lines) 1.5 M sequences. The contours display lev- 
els of constant J, labeled with log(//10 50 erg s). The values for 
the dotted contours are 3, 5, 7, and 9 times the value of the la- 
beled contour below them. From high to low, the six lines per 
sequence represent the following evolutionary stages: start main 
sequence (X c = 0.7), main sequence (X f ~ 0.25), start core He 
burning (Y c ~ 0.95), He burning (Y c ~ 0.4), start TP phase and 
during TP phase (>5 pulses). 
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Fig. 4. Integrated angular momentum J(M r ) = J Q j(m)dm di- 
vided by M^ 3 , as a function of the mass coordinate for the final 
stellar models of our sequences (cf. Table 1). The background 
contours display levels of constant J, as in Fig. 3. 



Ferrario & Wickramasinghe 2005) reaches up to the line con- 
necting 3 M Q and 15 M . The partly much slower rotation of the 
magnetic white dwarfs could be explained by additional mag- 
netic braking due to stellar winds in the white dwarf progenitors. 
However, our magnetic models still predict white dwarf rotation 
rates in the 1 M Q to 3 M range which are about one order of 
magnitude larger than what is found in pulsating white dwarfs 
(ZZ Ceti stars). This implies that the white dwarfs should lose 
additional angular momentum, either after they are born, or still 
inside their progenitor stars. 
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Fig. 5. Average core specific angular momentum of our initial 
and final models (cf. Table 1 and Fig. 4), versus initial stellar 
mass (full drawn lines). Points for 15 M stars computed with 
the same physics from Heger et al. (2005) have been added. The 
upper line corresponds to the initial models. Filled triangles cor- 
responds to the final models of the non-magnetic sequences, and 
filled squares to the final models of the magnetic sequences. The 
dashed horizontal line indicates the spectroscopic upper limit on 
the white dwarf spins obtained by Berger et al. (2005). Star sym- 
bols represent astroseismic measurements from ZZ Ceti stars 
(Bradley 1998, 2001; Dolez 2006; Handler 2001, Handler et al. 
2002, Kepler et al. 1995, Kleinmann et al. 1998, Winget et al. 
1994), where smaller symbols correspond to less certain mea- 
surements. The green hatched area is populated by magnetic 
white dwarfs (Ferrario & Wickramasinghe 2005; Brinkworth et 
al. 2007). The three black open pentagons correspond to the 
youngest Galactic neutron stars (Heger et al. 2005), while the 
green pentagon is thought to roughly correspond to magnetars 
(Camilo et al. 2007), where the vertical dotted green line indi- 
cates the possibility that magnetars are born with higher core 
angular momentum. 



In the first case, one may speculate whether the situation in 
analogous to the neutron star case: for pulsars it is well known 
that they spin down with time. As the ZZ Ceti stars have a cool- 
ing age of about 10 9 yr, they might have lost angular momen- 
tum on their cooling track. However, estimating the spin-down 
time through a magnetic wind by t\ - M^'-iM/BRXlGM/R^ 
(Justham et al. 2006) would allow for a significant effect at most 
in strongly magnetic white dwarfs, and spin-down through elec- 
tromagnetic radiation, with tt = Mc 3 l(B 2 QrR 4 ) is always in- 
significant. In the above expressions, M,R,Q and B are white 
dwarf mass, radius, surface magnetic field strength and spin fre- 
quency. 

In the second case, perhaps an additional internal angular 
momentum transport mechanism in white dwarf progenitors is 
required. While transport through gravity waves has been sug- 
gested as such (Zahn et al. 1997), we point out that, for stars 
above ~ 1.3 M Q , the core spin-down is likely to occur during 
the post-main sequence evolution, since on the main sequence 
these stars are still rapid rotators. It is unclear at present whether 
gravity waves can transport angular momentum through the shell 
source in the giant stage. On the other hand, the magnetic angu- 



lar momentum transport has been shown by Zahn et al. (2007) 
to be more complex than in the picture of Spruit (2002). Also, 
Yoon et al. (2008) find strong poloidal fields may be generated 
by convective cores, and suggest that their influence on angu- 
lar momentum transport may be comparable to that of the fields 
suggested by Spruit. 

5. Conclusions 

While nature may be more complex, it is worthwhile to attempt 
to understand the spins of compact stellar remnants with a sin- 
gle theoretical approach for all initial masses. We showed above 
that angular momentum transport through rotationally induced 
magnetic fields according to Spruit (2002) provides a major im- 
provement of the predicted spins of white dwarfs and neutron 
stars. Magnetic angular momentum transport appears also the 
most promissing candidate to bridge the still existing gap be- 
tween observed and predicted white dwarf spins at low mass (see 
Fig. 5). 

If internal magnetic torques in their progenitors are indeed 
resopnsible for the slow rotation of compact stars in the Milky 
Way, then there may be little room for an important role of angu- 
lar momentum in their formation process, at least in single stars. 
This may have implications for understanding the progenitors 
and the formation meachanism of magnetars (Sawai et al. 2008) 
and long gamma-ray bursts (Woosley and Bloom 2006). 
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